Variation of elastic scattering across a quantum well 
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The Drude scattering times of electrons in two subbands of a parabolic quantum well have been 
studied at constant electron sheet density and different positions of the electron distribution along 
the growth direction. The scattering times obtained by magnetotransport measurements decrease as 
the electrons are displaced towards the well edges, although the lowest-subband density increases. By 
comparing the measurements with calculations of the scattering times of a two-subband system, new 
information on the location of the relevant scatterers and the anisotropy of intersubband scattering 
is obtained. It is found that the scattering time of electrons in the lower subband depends sensitively 
on the position of the scatterers, which also explains the measured dependence of the scattering on 
the carrier density. The measurements indicate segregation of scatterers from the substrate side 
towards the quantum well during growth. 



The striking success of Ga[Al]As semiconductor het- 
erostructures originates from the extremely high mobil- 
ities obtained in these materials. One key ingredient 
for the fabrication of such samples is modulation dop- 
ing, where dopants and electrons are spatially separated. 
At low temperatures, impurity scattering, alloy scatter- 
ing and interface roughness scattering limit the electron 
mobility [Q. If more than one subband is occupied, in- 
tersubband scattering takes place in addition H[j| . 

Information on the relevant scattering processes is usu- 
ally obtained by measuring how quantum (r q ) and Drude 
scattering times (r) vary with carrier density ns. For 
two-dimensional electron gases (2DEGs) realized in Al- 
GaAs heterostructures, it is found that impurity scatter- 
ing is dominant. In this case, one finds r oc rig, with 7 
between 1 and 1.5, depending on the distance between 
the dopants and the 2DEG @]. 

In a two-subband system with subband densities ri\ 
and ri2, the Drude scattering times Ti of subband i are 
usually found to increase monotonically with [§,1). 
Recent results show that in a parabolic quantum well 
(PQW), Ti may also slowly decrease, i.e. 7 < 0, when a 
second subband is occupied [gj. In this paper, we inves- 
tigate this unusual dependence and show that it may be 
due to a certain arrangement of the ionized impurities. 

The PQW, grown by molecular beam epitaxy (MBE), 
is a 760 A wide Ah-Gai-^As layer with x varying parabol- 
ically between and 0.1 (inset of Fig. [l]a). In the cen- 
ter of the well, a three monolayer thick Alo.05Gao.g5As 
layer forms a potential spike. The well is embedded sym- 
metrically in 200 A of undoped Alo.3Gao.7As spacer layers 
and remote Si-doping layers on both sides. On the sur- 
face side, the donors are provided by 11 sheets, each with 
a Si donor density of nominally 5-10 m~ 2 
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FIG. 1. (a) Fit oip xx {B) for V /9 =-50mV (Vb g =+1000mV) 
to the two-subband scattering model. Inset: scheme of sam- 
ple layout along the growth direction, (b) Measured p xx 
for different electron positions along the growth direction at 
ns = 2.9 ■ 10 15 m -2 . Values for Ut g are indicated, and Ubg is 
varied between -2.2 V (top) and +2.2 V (bottom) in steps of 
0.4 V. Subsequent data are offset for clarity by 50 Q. From top 
to bottom, the electron distribution is displaced towards the 
substrate. The data for p xy fall on top of each other since ns 
is constant. Minima corresponding to the same filling factor 
in the lower subband are connected by a dashed line. 
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Si-concentration, arranged in a 200 A thick layer. On 
the substrate side, the donors are located within one 6- 
doping layer with a concentration of 5T0 15 m~ 2 . This 
asymmetry in the doping allows for saturation of the 
surface states and an effectively symmetric location of 
the electron distribution in the well. A back gate elec- 
trode consists of a 250 A thick n + -doped layer located 
1.35 (im below the well. A TiPtAu front gate electrode 
was evaporated on top of the structure. The experiments 
were carried out with standard Hall-bar geometries at 
temperatures of 100 mK. A magnetic field B was applied 
perpendicular to the electron gas. 

Figure |l|a shows a measurement of the magnetoresis- 
tivity p xx (B) at ns = 2.9 • 10 15 m~ 2 . From the low-field 
magnetoresistivity, T\ and Ti are obtained. Early studies 
on scattering times in two-subband systems relied on the 
assumption of two independent electronic systems with 
additive conductivities a = <j\ + a-x with Oi — riie 2 Ti/m, 
quantitatively explaining a measured positive magnetore- 
sistance [§],|],§| (e,m electron charge and effective mass). 
In a more sophisticated model based on the Boltzmann 
equation j^], intersubband scattering is taken explicitly 
into account. This leads to independent scattering times 



n(B) = Re J^(K + iud^kj/ki 



(1) 



where the ki are the Fermi wave vectors, fej = y27mi, 
c = eB/m, and K the scattering matrix defined by 
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The coefficients Pnm are related to the transition rates 
Pnmi'f) between subband states n and m and scattering 
angle <f> by Fourier transformation in <j). P^ is the tran- 
sition rate integrated over the allowed scattering vectors, 
while in pQ' the integrand is multiplied by cos <p. The 

difference P^ '' — P£' corresponds to the single-subband 
Drude scattering rate, where the matrix element of the 
scattering potential is weighted by (1 — cos#). Note that 
in the diagonal elements, also the isotropic part of inter- 
subband scattering is included. We have shown in a re- 
cent paper that intersubband scattering cannot be ne- 
glected in our experiments Q . 

With rii known, Eq. [l] allows a fit to p xx (B), with 
Ki,K 2 and K 3 being the fit parameters (||| (Fig. |l]a). 

We measured p xx (B) at ns = 2.9 • 10 m -2 (con- 
trolled by the low-field Hall voltage) and different po- 
sitions of the electron distribution along the growth di- 
rection (Fig. ^jp). The electrons were displaced by ap- 
plying voltages Uf g (t/b g ) between the front (back) gate 
electrode and the electron gas. 



Clearly visible are variations of both amplitude and 
period of the Shubnikov-de Haas (SdH) oscillations with 
changing Vf g . The amplitudes at a fixed magnetic field 
decay as the wave functions are displaced towards the 
substrate. This corresponds to a decreasing r q P,|lO||, an 
effect not to be discussed in this paper (see Rcf. [g for 
evaluated data on r q ). 

By fitting Uh s as a function of Uf g at constant ns to a 
capacitor model, we find the displacement Az per front 
gate voltage to be about 1000 A/V [0. Thus we can 
plot the data as a function of Az instead of gate voltages. 
From the SdH frequency we evaluate n\(Az) (Fig. ||a). A 
minimum occurs in n\ at U{ g « — 130 mV and is related 
to the narrow potential spike in the center of the PQW. 
The spike leads to subband energy shifts depending sen- 
sitively on the electron distribution along the growth di- 
rection. A displacement of the electrons thus changes ri\ 
and ri2. The difference between the two lowest subband 
energies reaches a minimum when the wave functions are 
centered with respect to the spike. Therefore, the mini- 
mum in ni provides the reference for the location of the 
wave functions in growth direction p2[ , where Az = 0. 

From the data, we evaluated t\ and t 2 for different Az 
(Fig. ||a) . Both t\ and T2 show a maximum as a function 
of Az. The maximum in t 2 occurs where wave functions 
are centered, i. e. Az = 0. 

Assuming a decrease of Tj with decreasing rii (7 > 0), 
we expect a minimum in t\ at Az = 0, which disagrees 
with the measurement. On the other hand, the scattering 
rate depends on the distances from the relevant scatterers 
p3| . For Az ~ 0, these distances are maximized, giving 
rise to large r^. The fact that t\ is large around Az = 
indicates that not its density-dependence dominates n, 
but the distance to the relevant scatterers. In contrast 
to the first subband, both, n,2 and T2 have a maximum at 
Az = 0. Note that n,2 is much smaller than m, leading to 
small Fermi wave numbers where screening is more effi- 
cient. Thus the screened scattering potential at relevant 
wave numbers is less sensitive to displacements along the 
growth direction. On the other hand, the relative change 
of ri2 with Az is larger than that of n\. Hence, T2 is 
stronger influenced by its density dependence than by 
Az, which explains the coincidence of the maximum in 
T2 with Az = 0. 

The maximum of t\ is shifted towards the surface, 
indicating stronger scattering on the substrate side. 
Although this could be explained by assuming more 
dopants than expected from the MBE growth protocol, 
we can excluded this, because the total amount of Si 
brought on the wafer was measured accurately. However 
there might be segregation of dopants on the substrate 
side towards the PQW during growth, which enhances 
scattering significantly. As we will show, a calculation of 
the n supports the assumption of segregated Si atoms. 
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FIG. 2. (a) Measurement of n, T2 and rti vs Az at 
ns = 2.9 • 10 15 m~ 2 . (b) Calculated scattering times with 
1.5 ■ 10 15 m~ 2 of dopants shifted towards the substrate-sided 
edge of the parabolic profile, (c) The same calculations as 
in (b), but with a distribution of scatterers as in the growth 
protocol. 



The matrix elements of the scattering potential were 
obtained by numerical integration using self-consistently 
calculated wave functions |14[| . Then the transition rates 

(i) 

Pnm were calculated by integrating the squared matrix 
elements over the allowed scattering vectors. Screening 
was included in the Thomas Fermi approximation. The 
Tj were calculated from Eq. [j]. A detailed calculation of 
the scattering rates based on different scattering mech- 
anisms reveals that the contributions of alloy scattering 
(including the potential spike) and interface roughness 
scattering are an order of magnitude smaller than that 
of Coulomb scattering. Initially, two layers of Coulomb 
scatteres were included. The dopants on the surface 
side were gathered in a single c5-layer 300 A above the 
well, with a concentration of N\ — 3 • 10 16 m~ 2 . The 
second layer is the doping layer 200 A below the well 
(N2 = 2.8 • 10 15 m -2 ). These values correspond to half of 
the nominal Si concentration brought on the wafer dur- 
ing the MBE-growth, qualitatively accounting for deep 
donors and not ionized impurities. Figure |^b shows the 
obtained scattering times. As expected for this donor 



configuration, t\ monotonically increases as the electrons 
are displaced towards the substrate side. 

In order to take segregated Si atoms into account, we 
placed A3 = 1.5 • 10 15 m~ 2 scatterers at the edge of the 
well on the substrate side, and reduced by the same 
amount (Fig ^c). As in the experiment, we obtain a 
maximum in t\ displaced towards the surface side and a 
maximum of T2 at Az = 0. At the surface side, t\ de- 
creases only slowly, saturating at a value comparable to 
the simulation with A3 = 0. It is the balance between 
the monotonically decrasing t\ shown in Fig. ^c, and the 
range and strength of the extra layer, which determines 
the exact shape of t\{Az) 
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FIG. 3. Measured (a) and catauftted (b) Ki,K 2 and -K 3 . 
In (b), ionized impurity scattering was modelled as in Fig.^c. 

The calculated scattering times are about 50 percent 
larger than the measured ones. It is well-known that 
for PQWs calculations overestimate the scattering times. 
Possible explanations are size-effect scattering from the 
edges of the electron gas or enhanced background im- 
purities due to the greater reactivity of Al with oxygen 
and carbon-containing molecules in the MBE chamber. 
In addition, the calculated values depend on how screen- 
ing of the scattering potential is implemented and which 
concentration of ionized impurities is assumed. We did 
not attempt to simulate accurately, here only the qual- 
itative behavior, in particular its spatial dependence, is 
of importance. 

Additional insight can be gained by studying the spa- 
tial variation of the matrix elements Ki (Fig. ||a). Usu- 
ally, Drude times are insensitive to small-angle scatter- 
ing. For intersubband scattering, K3 contains the part 
of the scattering rate weighted by cos</>. This gives in- 
formation about the amount of small-angle intersubband 
scattering. Since almost no structure in A3 is observed, 
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while K\ increases stronger on the substrate side, large- 
angle scattering must be higher on the substrate side. In 
order to increase large-angle scattering of Coulomb scat- 
terers with fixed density, the distance to the electron gas 
has to be diminished. This happens if scatterers segre- 
gate towards the electron gas. The calculated Ki nicely 
reproduce the experimental data jTig.||b). 





20 P 



10 - 



. (a) 




o° ° o oo oo- 




o 




c 


o- 

> 


V 

V 


_ 


o \ 

V T 2 











3 



(b) 



o 



v -V 







1 2 3 

FIG. 4. Measurements (^(ife'd Calculation (b) of scatter- 
ing times (symbols) and subband densities (lines) vs ns- In 
(b), small symbols are calculated without, large symbols with 
additional impurities at the substrate side of the well. 

With this strong evidence for segregated scatterers at 
the substrate side of the well, we come back to the previ- 
ously unexplained structure in the density-dependence of 
n §. In this experiment, Ub s was kept fixed, while Uf g 
and therefore n s was changed. In Fig. |I[ the measured 
and calculated values for 7i,T2,7ii and ni are shown. In 
the measurement, n slightly decreases as n,2 gets popu- 
lated. In the calculation, the additional scattering layer 
gives rise to a weak increase of T\ with ns when the sec- 
ond subband is occupied (large symbols), whereas a steep 
decrease results in the case of no additional layer (small 
symbols). Thus the additional scatters are responsible for 
the slope of Ti(ns). Since ns is driven by Vf g , the elec- 
tron distribution expands towards the surface side with 
increasing ns. Thus the scatterers on both sides of the 
well compete and determine the shape of T(ng). As dis- 
cussed above, for small n^, ti is not so much sensitive to 
additional scatterers, which is reflected in similar values 
obtained from the two simulations shown in Fig. ^b. 

In conclusion, we have presented an investigation of 
Drude scattering times in a modulation-doped multi- 
subband quantum well. Using front- and back gate volt- 



ages, the position of the electron distribution and the 
subband densities were tuned. The Drude scattering 
times of individual subbands were measured. It was 
found that t\ is dominated by the distance of the 2DEG 
to the impurities and not by its density dependence. Its 
behavior could therefore be used to locate additional 
scatterers at the substrate edge of the well, which are 
presumably due to segregation of dopants during growth. 
The measured scattering times could be qualitatively 
reproduced in a calculation assuming that half of the 
substrate-sided donors had diffused to the edge of the 
well. Using these results, previous measurements of the 
density dependence of t\ could be explained. While ob- 
tained for a PQW, the presented method of investigating 
the scattering times as a function of the electron-gas po- 
sition might give further informations on scattererers in 
other types of samples. 
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